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ABSTRACT 

We show that dissipationless ACDM simulations predict that the majority of the 
most massive subhaloes of the Milky Way are too dense to host any of its bright 
satellites (Ly > 10 5 Lq). These dark subhaloes have peak circular velocities at infall of 
VSnfaii = 30-70 kms -1 and infall masses of [0.2-4] x 10 10 M . Unless the Milky Way is 
a statistical anomaly, this implies that galaxy formation becomes effectively stochastic 
at these masses. This is in marked contrast to the well-established monotonic relation 
between galaxy luminosity and halo circular velocity (or halo mass) for more massive 
haloes. We show that at least two (and typically four) of these massive dark subhaloes 
are expected to produce a larger dark matter annihilation flux than Draco. It may 
be possible to circumvent these conclusions if baryonic feedback in dwarf satellites or 
different dark matter physics can reduce the central densities of massive subhaloes by 
order unity on a scale of 0.3 - 1 kpc. 
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1 INTRODUCTION 

The cold dark matter (CDM) paradigm has been demon- 
strably successful at explaining a variety of observations on 
cosmological scales. Tests on smaller scales are often compli- 
cated by the physics of galaxy formation, but are crucial for 
verifying the CDM model. Perhaps the most prominent issue 
facing ACDM on galactic scales is the large discrepancy be- 
tween the number of observed and expected satellite galax- 
ies of the Milky W ay (|Kauffmann et aL|l993| |Klypin et aL 



[19991 |Moore et aLl|1999| see |Bullock||2dlO[ for a recent re- 
view) . Accordingly, much theoretical work has been devoted 
to understanding how to reproduce the satellite population 
of the Milky Way (MW). 

While some ACDM models of the MW's satellite pop- 
ulation place the most luminous dwarf galaxies in the most 
massive subhaloes at redshift zero (|Stoehr et al.| |2002| 



Hayashi et al. 2003 Pcnarrubia et al. 2008 1, recent kinematic 



studies of the satellites have shown that this is unlikely to 
be the case ( |Walker et al.|2009| [Strigari et alT]|2010[ ). Other 
models postulate that MW satellite galaxies correspond to 
subhaloes that were the most massive at some earlier time 



(Bullock ct al. 2000 


Kravtsov et al.|2004b Ricotti & Gnedin 


2005 


Koposov et a 


.||2009| |Okamoto & Frenk||2009| |Busha 


et al.| 


2010 


), often the epoch of reionization, with galaxy 



formation strongly suppressed in lower mass subhaloes (see 
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Kravtsov|2010| for a recent review). In addition, many faint 
MW satellites have been discovered in the SDSS (e.g., Will- 
man et al~||2005 Belokurov et al.| |2007), and it has become 
clear that up to a factor of ~ 5 — 20 times as many faint 
galaxies could remain undetected at present owing to incom- 
plete sky coverage, luminosity bias, and surface brightness 
limits ( |Tollerud et al.|2008| [Walsh et al.|2009)|Bullock et al.| 
20101. 



While this theoretical and observational progress has 
alleviated - though not eliminated - concerns about the mis- 
match between the number of low-mass subhaloes and faint 
MW satellites, a pressing question remains: what is the stel- 
lar content of massive Milky Way subhaloes at redshift zero? 
In this Letter, we focus on properties of massive subhaloes in 
ACDM galaxy-mass haloes, and examine which Milky Way 
satellites - if any - can be hosted by such subhaloes. 



2 SIMULATIONS AND DATA 

Our ACDM predictions are based on dark matter subhaloes 



from the Aquarius project ( Springel et al. 2008 1 and the 
Via Lactea II simulation piemand et al.||2008[ |2007[ ). The 
Aquarius project consists of six galaxy-mass haloes - de- 
noted A through F - simulated at a series of increasingly 
high mass and force resolution. Although only halo A was 
simulated at the highest resolution level, all six haloes were 
simulated with particle mass m p = (0.64-1.4) x 10 4 M Q and 
Plummer-equivalent gravitational softening e = 66 pc; it is 
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this set of "level 2" simulations that we use in this paper. 
The Via Lactea II simulation (VL-II) of one galaxy-mass 
halo used m v = 4.1 x 10 3 M Q and e = 40 pc. A notable differ- 
ence between the simulations is the background cosmological 
model: the Aquarius simulations used a value of 0.9 for the 
power spectrum normalization erg and 1.0 for the spectral in- 
dex of the primordial power spectrum n s , while VL-II used 
as = 0.74 and n s = 0.951. The best current estimates of 
these parameters, cr 8 = 0.816±0.024 and n s = 0.968±0.012 
(Komatsu et al. 20111, fall in between those used for the 



simulations. 

In each simulation, we select every subhalo that lies 
within 300 kpc of the host's center and has a maximum 
circular velocity Vmax = max{[GA/(< R)/R] 1 ^ 2 } exceeding 
lOkrns -1 . We characterize a subhalo prior to infall onto its 
host via Vinfaii , which we define to be the value of V ma x when 
the subhalo's mass was at a maximum (over its entire evolu- 
tion) in Aquarius and the maximum value of V ma x over the 
subhalo's entire history for VL-II. The measured values of 
V m ax and i? m ax (the radius at which Vmax is attained) at red- 
shift zero are used to determine each subhalo's inner mass 
distribution by assuming that the subhalo's density struc- 



ture can be modeled by a Navarro, Frenk, & White ( 1997 



hereafter NFW) profile with the same V ma x and -Rmax. Us- 
ing subhaloes extracted from the Millennium-II Simulation 
( Boylan-Kolchin et al.]|2009[ |, we have verified that this ap- 
proach gives the correct mass to better than 10% at radii 
that are well resolvecf"| (as expected from earlier work by 
Kazantzidis et al.|2004 |. 

To connect the A^-body subhaloes to the bright (Lv > 
10 Lq) dwarf spheroidal galaxies of the Milky Way, we turn 
to kinematic measurements of the dwarfs' masses. IWalkerl 
|et al.| p009| and |Wolf et al.| |20Tof have recently shown 
that dispersion-supported galaxies such as the MW dwarf 
spheroidals have dynamical masses M1/2 within their de- 
projected half-light radii -R1/2 that are well-constrained by 
line-of-sight velocity measurements. Since these galaxies are 
all strongly dark matter-dominated even within -R1/2 (e.g., 
|Mateo|["l998[ ), observed values of M\/ 2 are effectively mea- 
surements of the dark matter mass within R±/2- A necessary, 
but not sufficient, condition for a subhalo to possibly host 
a given dwarf is that M sub (-Ri/2) agree with M dwar f (Ri/2). 
Conversely, a dwarf cannot live in a subhalo if A'f su b(-Ri/2) 
and Mdw a rf (R1/2) differ substantially. 
Given the values of M1/2 



calculated by 



(20101, we can therefore investigate what {-Rmax, V 



Wolf et al. 



ues of NFW subhaloes are consistent with the observed dy- 
namics of the bright MW dwarf spheroidals. We exclude 
Sagittarius, which is far from dynamical equilibrium, for 
the present. Figure [I] shows the resulting 1 a confidence re- 
gions in V mi x-flmax space for these nine dwarfs. The behav- 
ior of the contours for each of the dwarfs is qualitatively 
similar: there is a global minimum in V max, correspondin g 



to i? max = .R1/2 and V max = V^crios,* (Wolf et al 



2010) 



1 While both the host haloes and subhaloes from Aquarius are fit 
somewhat better by[Ei nasto ( 19 65|l profile s than by NFW profiles 
( |Navarro et al.|2010| |Springel et al.|2008^ , we use NFW profiles 
here because they provide more conservative constraints: at fixed 
V ma x and Rmax, an Einasto profile contains more mass than an 
NFW profile within a given radius R for reasonable values of the 
Einasto shape parameter a when R < i? max . 




Figure 1. Constraints on the V max — Rmax values (assuming 
NFW profiles) of the hosts of the nine bright (L v > 10 5 Lq) MW 
dwarf spheroidal galaxies. The colored bands show 1 a confidence 
intervals based on measured values of R1/2 and M]y 2 from IWolfl 
|et al.|p0Tu) . 



with allowed values of V ma x increasing for both smaller 
and larger values of R max (corresponding to Rmax < -R1/2 
and i?max > R1/2, respectively). The lower portions of the 
curves, where R max < -R1/2, are unlikely to be physically 
plausible models for the hosts of dwarfs, as they require that 
the dark matter subhalo has been strongly affected by tides 
on the scale of the luminous matter in the dwarf. 



3 RESULTS 

In Figure [2] we plot data for all subhaloes from the six 
Aquarius simulations (circles) and from the VL-II simulation 
(triangles), color-coded by Mnfa.11 • The gray shaded band cor- 
responds to 2 a constraints from the MW dwarf spheroidal 
galaxies in Fig. |T| In terms of the total mass within 300 par- 



sees (M300; Strigari et al.||2008| ), this gray shaded region is 
almost exactly the same as 6.5 x 10 6 < M 3O o/M < 3 x 10 7 . 
Many of the subhaloes lie in the range that is consistent at 
the 2cr level with the dwarfs, but there is a large population 
of subhaloes that does not. These subhaloes all have central 
densities that are too high to host any of the bright MW 
dwarf spheroidals; they also have higher values of both Vm ax 
and Vinf a ii, on average. 

The Milky Way contains three additional satellites 
that are brighter than the nine dwarf spheroidals included 
in Figs. [l][2] the Large Magellanic Cloud (LMC), the 



Small Magellanic Cloud (SMC), and the Sagittarius dwarf 
spheroidal. In the context of ACDM models of galaxy for- 
mation, the Magellanic Clouds are expected to reside in sub- 
haloes with large values of Vnfaii: using abundance matching 
(e.g., |Kravtsov et al.||2004a| |Conroy et aT1|2006| |Guo et al 
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Figure 2. Subhaloes from all six Aquarius simulations (circles) 
and Via Lactea II (triangles), color-coded according to Vi n f a ii . The 
shaded gray region shows the 2 a confidence interval for possible 
hosts of the bright MW dwarf spheroidals (see Fig.JTJ. 



2010 1 to assign stellar mass to subhaloe^J we find that the 
SMC should have VWaii = 70 - 80 kms -1 and the LMC 
should have Vinfaii = 95 — 105 kms -1 . Conservatively, we es- 
timate that the Magellanic Clou ds have Vnfaii > 60 k ms -1 
and Vmax(z = 0) > 40 kms -1 ( Stanimirovic et al 



2004 



Olsen & Massey 2007 \ and remove all subhaloes with these 



properties from our sample of subhaloes that are inconsis- 
tent with the dynamics of the bright MW dwarfs. Sagittarius 
is in the process of being completely disrupted by the disk, 
but estimates of its pre-interaction properties give it a to- 



tal stellar mass similar to the SMC ( Niederste-Ostholt et al 



20101. In the absence of the MW disk - e.g., in dissipation- 
less simulations such as those used here - it is very likely 
that Sagittarius would be much more massive at z = 0, and 
our Magellanic Cloud exclusion criteria might also be appro- 
priate for Sagittarius. Regardless, the inclusion or exclusion 
of one object does not alter the conclusions reached below. 
The remaining subhaloes are not compatible with hosting 
any of the bright (Ly > 10 5 Lq) satellites of the Milky Way; 
we refer to these as massive dark subhaloes and focus the 
remainder of our analysis on them. 

Figure [3] shows the cumulative velocity function of mas- 
sive dark subhaloes as a function of Vinfaii for each of the 
seven simulations considered here. All of the dark subhaloes 
plotted in Fig. [3] have current V ma x values larger than 23 
kms -1 , and none meet our criteria for hosting galaxies sim- 
ilar to the Magellanic Clouds. In all cases, there are at least 
6 - and up to 12 - subhaloes with Vi n f a n > 30 kms -1 that 
are not consistent with any of the bright MW satellites 
(i.e, any satellite with Ly > 10 5 Lq). These subhaloes tend 
to be more massive than the possible hosts of the dwarf 



2 We match n(> M*) from |Li fc White| pool?} to n(> V in fan) 
that we have calculated from the Millennium and Millcnnium-II 
simulations I Springel et al.|2005| [Boylan-Kolchin et al.|2009 1 . 



Figure 3. Cumulative Vinfaii function of massive subhaloes at z = 
that cannot host any MW satellite brighter than Ly = 10 5 Lq, 
including the Magellanic Clouds. Each of the seven high resolu- 
tion simulations studied here has at least six such subhaloes with 
Vnfall > 30 kms -1 , and at least four with Vi n f a li > 40 kms . 



spheroidals, both today and at infall (see Fig. [4] below). 
Moreover, the three haloes with the fewest massive dark sub- 
haloes (Aq-B, Aq-C, and VL-II) do not contain any potential 
Magellanic Cloud hosts. If we restrict ourselves to the simu- 
lations that do contain reasonable Magellanic Cloud analogs, 
then the predicted number of massive dark subhaloes is 
closer to 10, including several with Vinfaii > 50 kms -1 . 

The luminosity - Vinfaii relation for a representative halo 
is shown in Figure|4] The MW dwarfs (red symbols) were as- 
signed their Vinfaii values by placing the most luminous dwarf 
spheroidal (Fornax) in the subhalo with the largest value of 
Vinfaii that has M su b(-Ri/2) within 1 a of the measured Mi/ 2 
of Fornax, then repeating the process for each of the other 
dwarfs in order of decreasing luminosity. For each dwarf, 
the assigned value of Vnfaii can therefore be considered an 
upper limit at 68% confidence within this realization. The 
massive dark subhaloes (black symbols) are placed on the 
same plot according to their Vinfaii. These subhaloes must all 
have luminosities less than 10 5 Lq in order to have escaped 



detection in all-sky optical surveys (Whiting et al. 20071 



The dotted blue line shows an extrapolation of abundance 
matching, assuming M^jLy — 1. It is clear that neither the 
bright dwarf spheroidals nor the dark subhaloes described in 
this paper can be easily accommodated by galaxy formation 
models in which luminosity is a monotonic function of halo 
mass or V in faii- 

If massive, dark subhaloes do exist in the Milky Way 
halo, their presence has important implications for indirect 
dark matter searches. Denser subhaloes produce a larger lu- 
minosity from dark matter annihilation; from Fig. [2j the 
dark subhaloes expected in the Milky Way are denser than 
their potentially luminous counterparts and therefore may 
be bright in 7-rays due to annihilations. In Fig. [5] we plot 
the annihilation flux, J r , within 2.4 x 10~ 4 steradians (a cir- 
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Figure 4. Relation between \4 n f a u and Ly for Milky Way 
dwarf spheroidals (red points) and massive dark subhalocs (black 
points) for one representative halo realization (Aq-A). 



cular region with radius 0.5 degrees) of the center of each 
massive dark subhalo relative to the predicted flux within 
the same angular radius originating from Draco, one of the 
most promising targets among the MW dwarfs for Fermi 



(Abdo et al. 2010). The horizontal error bars on the data 



points show 68% confidence intervals based on 1000 ran- 
dom realizations for the observer's location (constrained to 
have a galactocentric distance of 8 kpc). Dark subhaloes 
are promising indirect detection candidates: each halo has 
at least two dark subhaloes with annihilation fluxes larger 
than that of Draco, and four of the seven haloes have at 
least one dark subhalo with J/Joraco > 5. 

Our division between dark and bright subhaloes is very 
conservative: rather than requiring at most one subhalo that 
is consistent with each of the bright MW dwarf spheroidals, 
we require that all of the dark subhaloes are inconsistent 
with all of the bright dwarfs at the 2 a level. While the 
quantitative results may, in principle, change slightly if sys- 
tematic errors in the determination Afj/a for the densest 
MW dwarfs (Draco and Ursa Minor) have resulted in an 
underestimate of M1/2, our general result - that many of 
the most massive simulated subhaloes are too dense to host 
any bright MW satellite - will be unaffected unless all of the 
measured M1/2 values change substantially. 



4 DISCUSSION 

The results of Section [3] show that high resolution ACDM 
simulations of individual galactic haloes generically predict 
that the Milky Way should host at least six subhaloes that, 
at one point, had maximum circular velocities in excess of 30 
kms -1 and yet are incompatible with any known MW satel- 
lite (including the Magellanic Clouds) having Ly > 10 s Lq . 
Either these subhaloes actually exist as predicted in the 
Milky Way, requiring us to understand their properties and 



Figure 5. Distribution of annihilation fluxes from dark sub- 
haloes, normalized to a typical scenario for the annihilation flux 
from Draco. Error bars reflect 68% confidence levels for varying 
the specific angular location of the observer on the solar circle. 
The typical halo has approximately four dark subhaloes with an- 
nihilation fluxes exceeding that of Draco. 



stellar content, or they do not exist, in which case we must 
understand the mechanism(s) that suppress their formation 
or survival. 

If massive dark subhaloes exist as predicted: 

Detecting massive dark subhaloes would be a strong 
confirmation of the standard ACDM paradigm. These dark 
subhaloes might host at least some of the recently discov- 
ered ultra-faint galaxies, all of which have luminosities lower 
than 10 Lq . Kinematic constraints favor masses and densi- 
ties for the ultra-faints that are indicative of fairly massive 



subhaloes (Strigari et al. 2008 Walker et al. 2009 Simon 



et al.|2uT0 l, albeit with large uncertainties at present (e.g. 
Wolf et al.| 2010||Martinez et al.|2010[ ). If some of the ultra- 



faints are hosted by the massive subhaloes described here, 
they would have total mass-to-light ratios of 10 5 — 10 8 . The 
ultra-faints would be excellent candidates for indirect dark 
matter detection in this scenario (Fig. [5]) . An alternate de- 
tection method could be through the subhaloes' tidal influ- 
ence on the MW's HI disk ( |Chakrabarti et al.|20lT| . While 
the existence of effectively dark subhaloes with low masses 
is perhaps not surprising given the standard ACDM power 
spectrum and the variety of effects that can impede cooling 
and star formation in shallow gravitational potential wells, 
the prospect of subhaloes more massive than the hosts of 
bright dwarf spheroidals but with Lv < 10 5 Lq is intrigu- 
ing. 

The existence of massive dark subhaloes requires that 
the fundamental assumption of abundance matching mod- 
els - that galaxy stellar mass or luminosity is a monotonic 
function of VWall _ does not hold for VWali < 50 kms -1 ; 
Fig. [4] illustrates this point. Galaxy formation on scales be- 
low 50 km s -1 should therefore be effectively stochastic, with 
stellar mass depending sensitively on specific details of a 
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subhalo's environment, formation history, etc. rather than 

primarily determined by host halo mass or VWali- 

If massive dark subhaloes do not exist as predicted: 

The most prosaic explanation is that the haloes studied 
here are not representative of the MW-mass halo popula- 
tion at large in ACDM. This is unlikely to be the case for 
the Aquarius haloes, however, since they have substructure 
abundances typical of the full sample of over 2000 MW-mass 
haloes from the Millennium-II Simulation's (137Mpc) 3 vol- 
ume ( jBoylan-Kolchin et al.|[2010l >. The subhalo mass func- 
tion of the Milky Way could also be a statistical anomaly 
with respect to ACDM expectations, in the sense that mas- 
sive MW subhaloes are all less concentrated than expected, 
or that there are zero massive dark subhaloes in the MW 
when we expect at least six (Fig. [3]). Perhaps the best way 
to investigate this possibility is to obtain detailed kinematic 
measurements of M31's satellites: if the M31 satellite system 
does not require massive dark subhaloes in ACDM simula- 
tions, then the statistical anomaly explanation would gain 
more traction. 

If the MW's subhalo mass function is not aberrant, 
then understanding why there are no massive dark subhaloes 
would likely result in important insight into the physics gov- 
erning structure formation. One possible ACDM-based ex- 
planation is that the dark matter distribution in satellites 
of the MW is substantially less concentrated than current 
dissipationless simulations predict^Baryonic processes may 
affect the dark matter distribution on small scales by heat- 
ing it to larger radii, which would have the desired effect 
of lowering the dark matter density. For such a solution to 
work, it would have to substantially lower the dark matter 
density on scales of 0.3 - 1 kpc (corresponding to the depro- 
jected half-light radii of the bright dwarf spheroidals) while 
not strongly impacting the dark matter on smaller scales 
(< 100 pc, corresponding to the half-light radii of the ultra- 
faint dwarfs), as ultra-faints seem to have high central dark 
matter densities ( Simon et al.|2010 1. It is not clear that this 
could produce nearly identical average dark matter densities 
on scales of 300 pc in galaxies spanning a factor of 10 4 in 
luminosity ( Strigari et al.|[2008 1 . 

Gravitational shocks from encounters with the MW disk 



VL-II and Aquarius bracket current estimates of erg and n„ . 
Modifying the dark matter power spectrum on sub-galactic 
scales - for example, through Warm Dark Matter (WDM) 
with a characteristic scale of 40 to 50 kms -1 - would result 
in both fewer massive subhaloes (e.g., Zavala et al. 20091 



may also destroy some fraction of satellites ( D'Onghia et al 
|2010[ ). This mechanism works most efficiently at destroying 
low-mass subhaloes, however. Furthermore, it would not af- 
fect subhaloes that have larger pericenters or were accreted 
recently; many massive dark subhaloes in the simulations 
studied here fall into these two categories. 

If the Milky Way's dark matter subhalo population is 
typical of ACDM predictions, and baryonic physics has not 
strongly modified the internal structure or abundance of 
massive subhaloes, then the more drastic solution of modi- 
fying the underlying cosmology may be required in order to 
circumvent our primary conclusions that massive dark sub- 
haloes should exist and that galaxy formation on small scales 
is stochastic. Merely tweaking the cosmological parameters 
within the ACDM model is unlikely to have an influence, as 



3 Lower concentrations of dark matter are also favored by many 
observations of low-mass field galaxies (e,g., |Kuzio de Naray et al.| 
|2008| l, although these tend to be gas-rich, disk-dominated systems 
with higher luminosities than the bright MW dwarf spheroidals. 



and lower central densities in such subhaloes. Recent anal- 
yses of the Ly-a forest put fairly stringent constraints on 
the mass of WDM particles, however ( Boyarsky et al.|2009 1. 
Dark matter self-interactions would also reduce the central 
densities of subhaloes, and would additionally make them 
more vulnerable to tidal disruption. It is far from obvious 
that the abundance and dynamics of observed MW satel- 
lites would be correctly reproduced in the viable parameter 
space of these non-CDM models. 

In summary, we find that the majority of the most mas- 
sive subhaloes in dissipationless ACDM simulations are too 
dense to host any of the bright Milky Way satellites. It fol- 
lows that galaxy formation must be effectively stochastic in 
haloes with maximum circular velocities of V < 50 kms -1 . 
This conclusion may be circumvented if the Milky Way's 
subhalo population differs substantially from the average 
ACDM expectation, or if the abundance or structure of mas- 
sive subhaloes in the Milky Way is strongly affected by bary- 
onic processes or different dark matter physics. 
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